Impurity-induced stabilization of solitons in arrays of 
parametrically driven nonlinear oscillators 



N.V. Alexeevalil, I.V. Barashenkovlll 

Max- Planck- Institut fiir Physik komplexer Systeme, Nothnitzer Str.38, Dresden, Germany 

G.P. TsironisS 

Physics Department, University of Crete and FORTH, P.O. Box 2208, 71003 Heraklion, Crete, Greece 

Chains of parametrically driven, damped pendula are known to support soliton-like clusters of 
in-phase motion which become unstable and seed spatiotemporal chaos for sufficiently large driving 
amplitudes. We show that the pinning of the soliton on a "long" impurity (a longer pendulum) 
expands dramatically its stability region whereas "short" defects simply repel solitons producing 
effective partition of the chain. We also show that defects may spontaneously nucleate solitons. 
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One of the most exciting developments in studies of 
collective synchronization of arrays of coupled nonlinear 
oscillators was the discovery of a dramatic increase of the 
degree of synchronization achieved by disordering the ar- 
ray . This raised hopes of the possibility of "taming 
spatiotemporal chaos with disorder" . In an attempt to 
understand the nature of this phenomenon, a numerical 
study of the effect of a single impurity on an otherwise 
homogeneous chain was carried out Surprisingly, it 
turned out that a single impurity can produce simple 
spatiotemporal patterns in place of complex chaotic be- 
haviour for very long chains of oscillators. 

The system simulated in Ref. Q (see also ) consisted 
of a chain of pendula coupled to their nearest neighbours 
and driven by a periodic torque. All individual pendula 
in the array were in their chaotic parameter regime while 
the natural frequency of the central pendulum (the im- 
purity) was out of the chaotic range. Consequently, the 
suggested mechanism of stabilization was the formation 
of a nonchaotic cluster around the defect which would 
subsequently pull the whole chain out of chaos (|,^. 

In the present note we study the effect of an impurity 
on a damped driven system by considering it from the 
viewpoint of nonlinear waves; that is, we explore collec- 
tive stabilization mechanisms. As in Ref. Q we study a 
chain of pendula but unlike [Q, we focus on the regime 
where all individual pendula are nonchaotic. (In spite of 
this fact the array as a whole may fall into the state of 
spatiotemporal chaos.) Another distinction from Ref. ^] 
is that our chain is driven parametrically not externally; 
the reason is the availability of explicit solutions. 

We assume that our pendula are strongly coupled. As 
a result, they tend to form soliton-like clusters of coher- 
ent behaviour. Similarly to externally driven systems, 
stationary solitons in parametrically driven chains are 
known to be stable for small driving strengths but lose 
their stability to oscillating solitons as the driver's ampli- 
tude is increased ||^-|^ . Increasing the driving amplitude 
still further, a spatiotemporal chaotic state sets in — 
with or without a series of intermediate bifurcations [|| . 



We will demonstrate that "long" impurities may act as 
centers of spontaneous nucleation of solitons and hence 
in chains with impurities solitons are even more generic 
structures than in homogeneous arrays. We will prove 
that pinning of a stationary soliton on a "long" impurity 
expands its region of stability. In particular, by choosing 
a sufficiently long impurity pendulum, the soliton can be 
stabilized in the parameter region where in the absence 
of the inhomogeneity it would set off the spatiotemporal 
chaos. Finally, although solitons pinned on "short" im- 
purities will turn out to be more exposed to oscillatory 
instabilities, we will show that such a pinning is an un- 
likely occurrence due to the repulsion between solitons 
and the short defects. 

The angle of the n-th pendulum in our chain satisfies 



— —mln {g + 4ijj^pcos2wt) sin0„. 



(1) 

where k is the torsion-spring constant, a the friction in 
the pivots, p and lo are the amplitude and frequency of 
the driver, and the length = 1 for all n ^ 0. In what 
follows we set g — m = 1. Defining the spacing between 
the pendula by a = l/Vk, the continuous limit of eq.(|l|) 
is given by the damped-driven sine-Gordon equation: 

[l+pSiz)]0rr-O,, + [l + qdiz)] X 

X [aOr + (1 + 4w V cos 2wt) sin 9] = 0. (2) 

Here 6{z) is Dirac's 6; p — «(?§ — 1) and q — a{lo — 1). 
Choosing the length of the central pendulum to be Iq — 
l±be with e <C 1, gives p = 2q = Aqe, where q = ±^ab. 

If the driving frequency is just below the edge of the 
continuous spectrum of linear waves: — 1 — e^, and as- 
suming that a and p are small: a — e^7, Auj'^p — 2he^, 
Eq.(||) has small-amplitude solutions of the form 

6 = 2etlj{T, X)e-"^^ + c.c. + 0{e^), 

where T = e^T 12, X = ez and V' satisfies the parametri- 
cally driven damped nonlinear Schrodinger equation 

iipT + Tpxx + 2|i/'| V - + 2q5{X)i; = h^p* - iji:. (3) 
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FIG. 1. A stable pinned soliton spontaneously created from 
a small- amplitude random spatial distribution. 

This nonlinear Schrodinger (NLS) equation was previ- 
ously used to model the nonlinear Faraday resonance in 
a long narrow water trough Q . The inhomogeneous term 
2q6{X)ip represents a local widening (for q > 0) or nar- 
rowing {q < 0) of the trough. The same equation de- 
scribes an easy-plane ferromagnet with a combination of 
a static and hf field in the easy plane and a biax- 
ial ferromagnet; in both cases the inhomogeneous term 
accounts for an impurity spin. Equation (j^) was also 
used in studies of the effect of phase-sensitive parametric 
amplifiers on solitons in optical fibers [l0|] . 

As in the spatially homogeneous case (q = 0) 
for h > yl + T^ ^^^"^ solution -0 = is unstable 
against continuous spectrum excitations. "Long" impu- 
rities {q > 0) harbour a discrete mode (here e <C 1): 



(5-0 = e cos nr 



l + h- 



-^sinr!T]e-^^-«l^l 



f]2 = (1 _ q2^2 _ ^i^g^ > ^^^^ ^ ^(X^^IjTq:^^ 
this localized mode also produces instability. For posi 
five q < -\/2, in which case f)^,^ 



lies below y/l + 7^ , the 



nonlinear development of this instability leads to the for- 
mation of solitons (figj^), stable or unstable. 

Equation (^ exhibits two stationary soliton solutions, 
each having a cusp at the origin: 



ip±{X) = ^±sech(A±|X| 



X) e 



-i9± 



where cos26i± = ±y^l-j^/h^, = 1 ± y/h^--/^, 
X = arctanh {q/ A±). For weak impurities, \q\ < 1, the ip+ 
soliton exists for any h and 7 satisfying /i > 7, whereas 
the requires, in addition, that h < l}q,-y. Strong im- 
purities {\q\ > 1) do not support the ip_ soliton at all 
whereas the 0+ exists only if h > i)q.'y. 

First we demonstrate that the soliton ijj- is unstable 
and can always be disregarded. Taking the linear pertur- 
bation in the form 5ip± — A±{f + ig)e~'^^^~'''^ gives 

- 5t - = Li/, ft - Tf = Lag, (5) 

where F = j/A^, t = A\T, and the operators 



-a^ + 1 - 6sech^(|a;| + i) 
1 =F 2i/ - 2sech^(|a;| + i) 



2QS{x), 
2Q5{x), 



(6) 
(7) 



72/Aj_, and 



with Q = q/A± = qVT^H, H = y^-'j-/^^, 
X = A±X. The minimum eigenvalue of the operator Lq 
associated with a nodeless eigenfunction sech(|a;| + x), is 
1^0 = ^2H. Consequently, in the case of the ip- soliton 
Lq is positive definite and Eqs.(m can be rewritten as 



(8) 



The maximum exponential growth rate A of solutions to 
EqM is given by In) 



= F^ 



</(.t)|-£i|/(.t)> 
7 < fix)\Lo'\f{x) > ■ 



sup 



(9) 



For any Q the operator Li has a negative eigenvalue 
^0 — 1 — associated with an even eigenfunction 



2/0 



(x) = e~'"^(3tanh^C + 3/ttanhf + - I), (10) 



where ^ = |a;| -I- 5; and k > 1 is a root of + 2k^Q + 
k{3Q^ - 5(5 - 4) -)- 3(5^ ^ 0. Thus the supremum in @ is 
positive, A is > F and the soliton 0_ is unstable against 
a symmetric nonoscillatory mode for all q, h and 7. 

A similar argument can be used to detect the disen- 
gagement instability of the other soliton, tjj^, arising for 
"short" impurities {q < 0). We simply notice that the 
second lowest eigenvalue of Lq which is associated with 
an odd eigenfunction 



wi{x) = sgn(2;)e'3(l^l+^) {tanh(|a;| 



Q}, 



is equal to fi = 1 ~ 2H - Q'^. For 2H < 1 - Q'^ or 
equivalently for h < i)q.-y, the operator Lq is positive def- 
inite on the subspace of odd functions. On the other 
hand, for Q < the operator Li has a negative eigen- 



value fii — 1 



associated with an odd eigenfunc- 



tion yi{x) — sgnx ya{x), with uq as in ( [ic| ) and k > 1 



a root of 



3kQ + 3Q —1 = 0. Hence the varia- 



tional principle (||) is still applicable, A is > F and the 
-0+ is unstable. The interpretation of this instability is 
straightforward if one notices that in the conservative 



(4) case (7 = 0) the inhomogeneous term —2qS{X)\^\^ pro- 
duces a local decrease respectively increase of energy for 
g > respectively q < 0. Consequently, in the conser- 
vative and weakly dissipative cases, the "long" impurity 
will attract and the "short" one repel small-amplitude 
tails of distant solitons. On the other hand, the energy 
of the pinned soliton is Eq = ^A^{1 - 3Q -I- 2Q^). For 
Q > (Q < 0) this is smaller (greater) than the en- 
ergy Eq of the infinitely remote soliton. These two facts 
indicate that Q > 0-impurities should attract and trap 
solitons (cf. [^). In the Q < case, conversely, distant 
solitons will be repelled while an initially pinned soliton 
will depin and move away from the impurity regaining 
its cusp-free shape. 
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FIG. 2. The effect of the g < impurity. For h < 
the growth of the asymmetric instabihty leads to the soU- 
ton's depinning and repulsion from the impurity (a), while 
for h > i)q,-y, the asymmetric Hopf instability evolves into a 
two- humped soliton staggering in a seesaw fashion (b). 

This was indeed confirmed by simulations of Eq.(|^) 
(fig.|^(a)). It is fitting to note that the depinning in- 
stabihty is not connected with overdriving the chain; it 
occurs already in the undriven NLS . 

In the region h > l)q^j (as well as in the case of sym- 
metric instabilities, and for long impurities q > 0) the 
variational principle ^ is not applicable. Here we let 
f{x,t) = w(x)e*^* and g{x,t) — — p^-^ f (x)e'"*, where 
= ^2 _ p2^ reduces to an eigenvalue problem 

LiU — UV, LqV — LOU, (11) 

which we solved numerically. Notice that we have re- 
duced a three- to two-parametric problem. Having found 
Lu = Lu{q,H), one immediately recovers the instability 
growth rate |Imf2(g, H)\ — T for all q, h, and 7. 

For small h < i)q^^ there is only one unstable pair of 
imaginary eigenvalues ±uJi, with the associated u and v 
being odd. As h is increased beyond ijq,^, the imaginary 
eigenvalues move onto the real axis (i.e. the soliton resta- 
bilizes) while another real doublet ±1^2 detaches from the 
continuous spectrum. The two collide and emerge as a 
complex quadruplet after which the real and imaginary 
parts grow until an asymmetic instability sets in. This is 
now not a disengagement instability; the stationary soli- 
ton is replaced by a two-humped structure (still pinned 
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FIG. 3. Atlas of stability charts on the (7, /i)-plane. Above 
the dashed "horizontal" curve, h — all localized so- 

lutions are unstable w.r.t. continuous spectrum waves. Solid 
"horizontal" curves are h = i)q.^. Below their corresponding 
flq, ^-curves impurities with q < repel solitons. The family 
of "parabolas" depict the onset of (symmetric) instability of 
the pinned soliton; the greater is q the larger is the stability 
domain. Finally, in the region between + and f)^,^ a 
q > 0-impurity will spontaneously nucleate solitons. 

on the impurity) whose left and right wings oscillate 180° 
out of phase (fig.^(b)). 

When q > 0, the motion of eigenvalues lu on the 
complex plane is similar to the homogeneous case [|]. 
The soliton -0+ is stable for h close to 7 but loses its 
stability to a symmetric oscillating soliton as h is in- 
creased. Fig.|| shows the Hopf bifurcation curves h = 
hg{'j) obtained from the relation \lmn{q, H)\ = T for 
q = 0.1, 0.3, 0.4, 0.5 and -0.3. For q > the stability 
domain is wider than without an impurity. For example, 
in a chain with the coupling k = 0.5 driven at lu ^ 0.98, 
lengthening the central pendulum by 15% (which gives 
q — 0.5) is sufficient to double the size of the soliton's sta- 
bility domain. If the same chain is driven at w ~ 0.999, 
the q = 0.5 and hence the same effect is produced just by 
a 3% elongation of the central pendulum. On the con- 
trary, the q < impurity narrows the stability domain. 

Thus, long impurities exert a strong infiucnce on soli- 
tons' dynamics. For h < t)q .y they attract and trap 
solitons; for h > t)q_-y pinned solitons are spontaneously 
formed around the q > defects. On the other hand, 
the soliton with h and 7 such that it would ignite spa- 
tiotemporal chaos in the homogeneous case is stabi- 
lized when pinned on a sufficiently long impurity (figj^). 
Therefore the q > defects should have a stabilizing ef- 
fect on the chain. One should keep in mind, however, 
that spatiotemporal chaotic states are not localized and 
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FIG. 4. In the homogeneous chain with 7 = 0.315, h = 0.95 
the unstable soliton seeds spatiotemporal chaos (a) . Introduc- 
ing the impurity with q = 0.3 is sufficient to stabilize it (b). 



a single stable soliton will clearly be insufficient to sup- 
press chaos in a long chain. The chaos can always be trig- 
gered by choosing the initial condition far enough from 
the soliton. In order to suppress chaos in a larger phase 
volume multiple impurities should be introduced; one is 
therefore led to the necessity of examining stability of 
solitons and periodic waves on finite intervals. 

Next, the fact that short impurities enhance the sym- 
metric instability should not play a destabilizing role 
since solitons tend to avoid "short" defects. On the 
contrary, these repulsive inhomogeneities will effectively 
partition the chain into smaller intervals and this will 
generally have a stabilizingcffect since long-wavelength 
instabilities will not fit in |]l3| . 

The analysis of the effect of multiple impurities (as well 
as the related case of finite intervals) is beyond the scope 
of this work. Here we only mention that introducing 
more than one impurity can result in more complicated 
(though still regular) patterns. We illustrate this by sim- 
ulating the case of two impurities. As one could expect, 
when two stable stationary solitons are pinned very far 
from each other, they do not interact and remain time- 
independent. However, if the separation is smaller than a 
certain critical distance, they start exchanging weak radi- 
ation waves and develope spontaneous oscillations which 
subsequently synchronize (fig.^.) 
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FIG. 5. Two q > impurities placed in the middle and at 
the end of the periodic interval, respectively. 
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